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Abstract 

The formalism for analyzing the magnetic field distribution in the vor- 
tex lattice of Pauli limited heavy electron superconductors is applied to 
the evaluation of the vortex lattice static linewidth relevant to Muon Spin 
Rotation (^iSR) experiment. On the basis of the Ginzburg-Landau ex- 
pansion for the superconductor free energy we study the evolution with 
respect to external field of the static linewidth both in the limit of indepen- 
dent vortices (low magnetic field) and in the near H^ 2 (T) regime by using 
an extension of the Abrikosov analysis to Pauli limited superconductors. 
We conclude that in the Ginzburg-Landau regime with Pauli limit, the 
electrodynamics of the vortex lattice predicts anomalous variations with 
applied field of the static linewidth which is a result of the spin response 
contribution to screening supercurrents that dominates the usual charge 
response. The model is proposed as a benchmark for comparison with 
possible other effects including vortex core localized states or interplay 
with magnetism. 

1 Introduction 

As an example of a superconductor in the Pauli limit, the heavy electron sys- 
tem CeCoIns has special properties regarding its response to external magnetic 
field. Notably Muon Spin Rotation /iSR experiment pQ has revealed anomalous 
variations of the vortex lattice static linewidth o\^ L with respect to magnetic 
field oriented along the tetragonal crystal c-axis. In Ref. pQ the static linewidth 
was measured at temperature T = 20mK, showed an increase with applied field 
from zero field to about 95% the upper critical field and eventually decreased 
just before the first order superconductor to metal transition. The decrease 
in <7^ L with respect to external field usually observed and analysed [2] is the 
hallmark of a diminution in the vortex lattice local field contrast due to the 
decrease of the inter-vortex spacing with increasing field. 

Here it is shown that in the Ginzburg-Landau regime an increasing behavior 
of the static linewidth is predicted. This results from the Zeeman interaction of 
the electron spin with the superconductor internal field, which dominates the 



usual charge response supercurrents. As a result the field distribution is modified 
on a distance ~ £ (£ is the coherence length obtained in the Ginzburg-Landau 
formulation) from the center of each vortex [3] . The existence of the effect was 
pointed out [I] in the context of magnetism of the FFLO (Fulde-Ferrel-Larkin- 
Ovchinnikov) state. In parallel to this, a numerical approach to Eilenberger 
equations was undertaken [5] and effects of strong Pauli paramagnetism were 
highlighted in the vortex lattice state of Pauli limited superconductors. 

Let us first discuss qualitatively the properties of Pauli limited heavy elec- 
tron superconductors [H [3]. This class is characterized by a zero tempera- 
ture Maki parameter here defined ajv/o = H° 2 q/H^ 20 > 1 ( an alternative def- 
inition includes a \/2 factor which is not assumed here for clarity). We set 
H°2Q = (f>o/£,o an d Hc20 = ^c/Mj tne zer0 temperature scales for orbital and 
Pauli limiting fields respectively (we use throughout units where K = c = 1). 
4>q = 7r/e ~ 2.07 x 10 7 G.cm 2 is the vortex fluxoid quantum, e the absolute value 
of the electron charge, £o = v f/T c the T = Cooper pair radius or coherence 
length, vp — kp/m* the Fermi velocity, kp the Fermi momentum, m* the renor- 
malized electron mass, T c the superconductor critical temperature, /i = g/is/2 
the electron magnetic moment absolute value, g the Lande factor, [Ib — e/(2m) 
the Bohr magneton, and m the electron bare mass. 

There are three characteristic lengths in the problem: the zero temperature 
coherence length £ defined above, L(H^ 20 ) = ^fi(f) /T c the inter-vortex dis- 
tance of a square vortex lattice in the Pauli limit at temperature T = and 
field H^ 20 (more generally we note L{B) = ^4> /B the inter-vortex spacing of a 
square vortex lattice with internal field B), and the London penetration depth 

= yj m*/(47rne 2 ) with n the electron density in the superconductor (at T=0 
and for a cylindrical fermi surface this is the electron density in the 2D metal 
n = k F /(2nl c ) with l c the spacing between the planes of the tetragonal crystal). 
Hence 



«M0 



^4>qT c m*T c 



and we define the Ginzburg-Landau ratio 
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(2) 



where r e = e 2 /m is the classical radius of the electron and r e kp ~ 10 -5 . 

The orders of magnitudes are as follows. In a classical, non-heavy electron 
superconductor, m* ~ m and Ep ~ I0 3 T c give k ~ 1 and omo ~ 10~ 3 . In 
CeCoIn 5 however, T c ~ 2.3K, £ - fOOA, X L ~ 50001 yield m* - 300m, 
Ep ~ I00T c , k ~ 50, and cxmo ~ 3, which is as we shall see the origin of 
special magnetic properties of the vortex lattice. The large Ginzburg-Landau 
parameter implies [5] at T = the ratio between the field at which the first 
vortex nucleates in the bulk of the sample and the orbital upper critical field 
H c w/H°2Q ~ \u(k)/k 2 <§C 1, hence B ~ H for a broad magnetic field range. In 
a Pauli limited superconductor, H cW /H^ 2f) ~ (r e kp)(Ep/T c )\n(K) ~ I0 -3 and 
the same property applies. 
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We now study the electrodynamics of the vortex lattice which results from 
such large values for parameters and ([2| . The vortex lattice static linewidth 
is defined as 

VL _ 7m 



a^ = ^5h(v)*, (3) 

where 7 M = 2tt x 135.5342MHz/T is the muon gyromagnetic ratio, h(r) is the 
component of the internal local field parallel to the applied field H, 6h(r) = 
h(r) — B, the macroscopic internal field (or induction) B — h(r), and overlinc 
means averaging over a vortex lattice unit cell. Eq. ^ can be expressed as a 
sum involving all order Fourier components F mn of the field distribution in the 
vortex lattice 

VL _ 7m 



E (^«) 2 - ( 4 ) 

(m,n)/(0,0) 



The components F mn are called vortex lattice form factors [3J |7j in the context 
of Small Angle Neutron Scattering (SANS) experiment [8]. 

2 Muon static linewidth in the low-fleld-high- 
temperature regime 

Here we use results of the Ginzburg-Landau formulation [U |3] to evaluate the 
static linewidth Eq. Q. The near-T c Ginzburg-Landau regime in the Pauli 
limit is accessible since the crossover temperature T* from orbitally limited to 
Pauli limited superconductivity is in the range (T c — T* ) /T c ~ 1 /a\ I0 (we shall 
see below this follows from the relation H°^° (T) / ff c p 2 (T) - a M0 ^fl - T/T c ). In 
the independent vortex approximation (low magnetic field) and high-K limit the 
form factors can be decomposed as a sum of two distinct contributions [3 

r mn — r mn • £ mn- \° I 

The first term is the usual charge response which gives rise to orbital supercur- 
rents. This writes [3] 

FZn = (6) 

where = v2£ is a variational parameter which minimizes the superconductor 
free energy q rnn = [2n / L(B)](m 2 -fn 2 ) 1 / 2 in a square vortex lattice, and K n (z) 
is the n th order modified Bessel function of the second kind (or MacDonald 
function) [5]. The near-T c coherence length and penetration depth depend on 
the symmetry of the superconducting gap [3] . The expressions for d-wave pairing 
are 



1 32tt 2 T 2 



e Kim 



and 

1 16 2 2 »r 

^2 = ~J nC VF 



Tr. " V 2ttT, 
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Figure 1: Variations of 



the dimensionless /iSR static linewidth 
■J L is taken from Eq. Q and the form 
factors from Eq. (|ll|). Different values for the temperature Maki parameter 



y2(27rA)7( 7 ^ ) where 



(Il3|) where used as indicated in the legend. 



where N is the density of states with dimension [Energy x Volume] -1 , £(z) is 
the Riemann zeta function, £(3) ~ 1.20. The field at which these two lengths 
diverge is defined as the near-T c Pauli limit upper critical field 



H? 2 = 2n'^J Tc , T . (9) 

c2 n y 7C(3)T C V ' 

The near-T c Zeeman spin contribution in Eq. |5]) for superconducting gap with 
d-wave symmetry reads [3] 

F mn = 3^ ( X ) K (q mn S v ). (10) 

We scale the internal field, form factors and coherence length such that |5]) in 
dimensionless units evaluates to 

fmn = 2 q 2 K 1 (q)+47rb 2 K Q {q). (11) 
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Here b = BjH% 



q = q mn i v = J^/ffl — ^ + ^ ( 12 ) 
V 2 a M 

with 

Q!M = ot M Q\J 1 - ^r, (13) 
«mo is given in Eq. (flj), and / mn = F mn {2-K\) 2 /<f>Q. Near T c the dimension- 



less form factors (11) take on a simple, universal form where only remains the 
parameter olm controlling the relative contributions of the spin response with 
respect to the charge response. The static linewidth variations in the indepen- 
dent vortex limit with dimensionless internal field b and different values for ajf 
are shown in Fig. [T] Observe the low field regime were all curves meet which 
follows from the limit 



^ ^ (/mn) 



2 -> 



Z_— / | ^ -U 71 Z I ^ 



V ^— ' \ ^— ' (m 2 +n 2 ) 2 

y (m,n)/(0,0) \ (m,n)/(0,0) v ; 

as 6 -> 0. (14) 

We now turn to the effect of temperature on the form factors Eq. ([5]). We fix 
the value a M o = 3 and plot <j^V2[2ir\{0)] 2 /f/y^o) for different T/T c , where 
A(0) is the Ginzburg-Landau penetration depth Eq. ^ taken at T = 0. The 
results are shown in Fig. [2j Notice we have extended the temperature domain 
to very low T/T c , which is not formally justified as it is in the near-T c region 
but is expected to give qualitatively meaningful variations. 

The MacDonald functions assume the limits K (q) — > —hx(q/2) — C and 
Ki(q) ^l/(|as(j->0 where C ~ 0.5772 is the Euler constant. Having SANS 
experiment in mind in the large-ckM limit it is useful considering Eq. ( |11[ ) with 
(m,n) = (1,0) 

W-M^SgStfo). (15, 

3 Muon static linewidth close to the second- 
order transition critical field 

The form factors as given in Eqs. ([§ and (jToJ are found in the independent 
vortex approximation. The derivation does not work near the transition to the 
non-superconducting metal. In the high-field limit close to the transition line, 
the main source of magnetic field inhomogeneity in the vortex lattice comes 
from the Zeeman spin response |H |3] 

5h(r) = -47re(\A(x,y)\ 2 -\M^yjf), (16) 
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T/T c =0 
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Figure 2: Variations of the static linewidth cr,^ L -\/2[27rA(0)] 2 /(7 M <?!>o) with field 
at different temperatures. We considered the parameter ccmo = 3. 



where 

by overlining we again mean averaging over a vortex lattice unit cell, and \|>( n ) (z) 
is the polygamma function [5] of order n. In a square vortex lattice the Fourier 
decomposition of the square of the gap magnitude reads [3] 

\A(x,y)\ 2 = \A(x,y)\* £ { _ ir+ n +mn 

m,n— — oo 

xe ~f (m 2 +n 2 ) e 2irimx/a e 2Triny/a ^g'j 



Therefore the form factors corresponding to Bragg peaks with indices (m, n) ^ 
(0, 0) take on the form 



F mn = -^e\A{x iy W{-l) m+n+mn e-^ m +" \ (19) 
and the vortex lattice static linewidth simply reads 

CTsVL = ^^|A(x,y)| 2 , (20) 
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T = 1.3K 




Figure 3: /iSR static linewidth close to the second-order transition line H^ 2 (T) as 
obtained from Abrikosov's analysis in the Pauli limit Eq. p0| ) for a temperature 
range as indicated in the legend. We have scaled the internal field with respect 
to H^ 2 (T) here defined as the curve solution of a(T,B) = [4] [3] . Notice the 
rapid increase in the absolute value of the slope of a^ h (B) while approaching 
the first order transition at T/T c ~ 0.56 and /J,H* 2 /T C ~ 1.07. 



where 



e- m{1 ) - 1 - 0.4247. 



(21) 



Eq. (20) shows explicitly that the vortex lattice contribution to the static 



linewidth vanishes when the transition is of the second order but shows a dis- 
continuity where the transition is of the first order. In the former case, the gap 
average is [H |3] 

W^W=^- a - (22) 



2/?a/3 



Here 



a = N 



In 



2 2ttT J 



(23) 
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and 



64tt 2 T 2 



37V C 



( 



2 2nT 



1 . \iB 



) 



(24) 



are the quadratic and quartic coefficients of the Ginzburg-Landau free energy [3J 

2 

[3] respectively, and ^l\z) is the digamma function 9J- Pa = y)| 4 /|A(a;, y)| 2 

is the Abrikosov parameter, it is = 1.18 for a square vortex lattice and 
/3 A = 1.16 for a triangular lattice. Then it follows 



4 Conclusion 

On the basis of the Ginzburg-Landau expansion for the superconductor free 
energy in the Pauli limit, we have studied the evolution with respect to ex- 
ternal field of the muon spin rotation vortex lattice static linewidth both in 
the limit of independent vortices (low magnetic field) near T c , and in the near 
H^ 2 (T) regime. In the first case, we have found a simple form for the total 
form factor which is a function of the internal field scaled with the temperature 
dependent upper critical field in the Pauli limit and includes a single parameter 
dim = a a/0 \A ~~ T/T c with a mo = ^4>qT c /v 2 f . In the regime near H^ 2 (T) we 
have applied an extension of the Abrikosov analysis to Pauli limited supercon- 
ductivity and observed the approach from second-order to first-order transition 
to the metal (this occurs at T/T c ~ 0.56K and fiH^ 2 /T c ~ 1.07) with a sudden 
raise in the absolute value of the slope of cr^ L (S) while approaching H^ 2 (T). 
Such an analysis allows a simple modelling of the effect of heavy electron su- 
perconductor strong paramagnetism on the vortex lattice electrodynamics. It 
is proposed as a benchmark for studying new puzzling vortex lattice properties 
in CeCoIn 5 [10]. 
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